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CALCULATION OF AERODYNAMIC LOADING AND TWIST CHARACTERISTICS 
OF A FLEXIBLE WING AT MACH NUMBERS APPROACHING 1.0 AND 
COMPARISON WITH EXPERIMENT 

By John P. Mugler, Jr. 


SUMMARY 

An iteration method is presented by which the 
detailed aerodynamic loading and twist characteris- 
tics of a flexible winy with kno wn elastic properties 
may be calculated. The method is applicable at 
Mach numbers approaching 1.0 as well as at sub- 
sonic Mach numbers. ( nictitations are made for 
a winy-body combination; the winy was swept back 
,yrJ° and had an aspect ratio of 4 . Comparisons 
are made with ex per' mental results at Mach numbers 
from 0.80 to 0.08. The calculated and experimental 
loadinys and twist distributions are in yood ayree- 
ment at a Mach number of 0.80 at low anyles of 
attack. With increases in anyle of attack or Mach 
number , separation and mixed-flow phenomena cause 
the a yree me nt between the calculated and experimental 
results to become poor . 

INTRODUCTION 

The aerodynamic loading on a flexible lifting 
surface is greatly affected by the structural 
deformations caused by the loading. Therefore, 
the calculation of the aerodynamic loading on a 
flexible surface is complicated by the fact that 
the elastic characteristics of the surface must be 
taken into consideration in order to obtain a good 
estimate of the resultant aerodynamic loading. 
References 1 to 3 present methods developed 
for the calculation of the aerodynamic loading 
on flexible wings. These methods are based on 
lifting-line concepts and a detailed theoretical 
description of the flow cannot be obtained. Since 
the boundary conditions for the lifting-line theory 


are satisfied at only one point on the chord, 
aeroelastic deformations in the chordwise direction 
cannot be taken into account in the calculations. 

Presented herein is an iterative method, based 
on lifting-surface concepts, for the calculation of 
the aerodynamic loading on flexible- wings. For 
lifting-surface theory, the boundary conditions are 
satisfied at several points on the chord, and the 
effects of chordwise deformations are accounted 
for in the calculations. The results yield both 
a detailed theoretical description of the flow on 
the flexible surface and the shape of the aero- 
elastically distorted surface. 

Calculations were made by the method de- 
scribed herein for a wing-body combination; the 
wing was tapered and had a low aspect ratio. 
Theoretical and experimental results arc com- 
pared at Mach numbers from 0.80 to 0.98. 

SYMBOLS 

a nm coefficients in expression for lift 

c section chord 

c wing mean aerodynamic chord 

c av average chord 

Ci section lift coefficient, — 

' qc 

Tfl 

c m section pitching-moment coefficient. — * 

1 & ’ qcr 

I) diameter 

K total down wash factor 

f dowmvash factors for trailing vortex 

G nv chordwise replacement load 

y dowmvash factors for bound vortex 

K kernel function 
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/ 

section lift, 1 A Idx 



lift at any point 


.1/ 

freo-stream Mach number 


m 

section pitching moment about 

c/4, 




A/> 

difference between lower-surface 

pres- 


sure and upper-surface pressure 


q freo-stream dynamic pressure 

/• body radius 

7 dimensionless body radius with respect 

to .Vp 

S wing area 

n wing semispan 

s' unsupported semispan (lateral distance 

from outer face of wing mounting 

block to tip) 



semispan of vortex 

t 

tangent of angle of sweep of load lines 

V 

freo-stream velocity 

w 

dow nwash velocity 

X 

distance from midchord line to lift 
point 

V 

longitudinal and lateral coordinates, 
respectively 

x, 7 

dimensionless coordinates, xjs t and y/$ v , 
n*speetively 

jjp 

wing-section center of pressure meas- 

c 

ured from leading edge 

V 9 

lateral distance from outer face of wing 
mounting block 

a 

angle of attack 

At* 

difference between angle of attack of 
wing station and angle of attack of 
w ing-body center line 


~aF 

V 

dimensionless coordinate, y/s 

v 9 

dimensionless coordinate, y'l* f 


0=cos 1 A: 
e/2 

C tj wing-twist influence coefficients due to 

normal loads at c/4 

D (J wing-twist influence coefficients due to 

moments about c/4 

Subscripts: 

a index of control points 

b index of lift point or vortex 


Matrix notation: 

{ } column matrix 

[ ] square matrix 

[ j _1 reciprocal matrix 

METHOD 

The method presented herein may be outlined 
as follows : 

(1) Calculate the aerodynamic loading on the 
rigid wing using lifting-surface theory. The 
known built-in angle-of-attack distribution over 
the sur ace is used to define the boundary condi- 
tions. 

(2) Use the elastic properties of the surface in 
conjunction with the aerodynamic loading from 
step (1 to calculate a new angle-of-attack distri- 
bution over the surface (corrected for aeroelastic 
deflect! >ns). 

(3) hep eat steps (1) (using the corrected angle- 
of-attack distribution) and (2) until the aerody- 
namic and structural forces are in equilibrium. 
The results yield both the detailed aerodynamic 
loading on the aeroelastieally distorted surface and 
the surface shape under load. 

AERODYNAMIC LOADING ON RIGID WING 

The ift ing-surface theory outlined in reference 
4 was used to calculate the aerodynamic loading 
on the rigid wing. The integral equation that re- 
lates a known downwash distribution to an un- 
known lift distribution may be written as 

7/?(.r U} y a ^ ^ // b ) X •Ffn 1/b }ja)d J'b* b 

Li equation (1), w(x a ,y fI ) is the downwash velocity 
at the )oint ; K is the kernel function and 

re prose its the contribution to the downwash at the 
point (c a ,y a ); and A l(x b ,y b ) is the unknow n lift dis- 
tribution or local doublet strength associated wil h 
a press ire doublet of unit strength located at any 
point (r ft5 7/ ft ). In reference 4 the wing surface is 
replaced by a network of horseshoe vortices. The 
unknown loading is assumed to be expressible as a 
series ( f loading terms each weighted bv an un- 
known coefficient to be determined. The series 
expression for the loading is used to obtain a set of 
chordwise replacement loads to represent the 
strength of each vortex. The downwash per unit. 
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load produced by each vortex is expressed analytically as a dowmvash factor. Through this proce- 
dure the integral equation can be approximated by the summation (see eq. (14) of ref. 4) 


Ac, \r 9 S '*, q r 

w(x a , y a )=^- S XI S Gnrfijl- 

7T( a= 1 j/=l //=<) 


- V 2 ( rt n(l+ <r/ »l r ?+ • • • a nm 7 i m )fi'(‘ r a)M>Qyf ]h V a) 


( 2 ) 


where a„ m represents the unknown coefficients, G nv , the replacement loads, and E, the dowmvash 
factors. In order to determine the values of the unknown coefficients a ntn , a set of control points is 
selected equal to the number of unknown coefficients. Then, if equation (2) is written for each 
control point and if the dowmvash angle wj\ T is equated to the slope of the mean camber line at 
each control point, a set of linear simultaneous equations results from which the unknown coeffi- 
cients &nm can be determined. This set of simultaneous equations can be written in matrix formas 


{n}=[K\{A} 


( 3 ) 


In equation (3) the elements of the column matrix [B] are the slopes of the mean camber line at each 
control point and are defined as 


tan ai 
tan 


{B} = 


tan olj\ 


(4) 


The elements of the square matrix [E] arc the products of the dowmvash factors and the replacement 
loads and are defined as 


w 


4*r 

7 rC 


£ 1Z Q*+f\ -tfEfa.MAfj'-yi) £ £ ,MA,V*-yi)n 


Q 


#/• . Q 


9,(9 r Q 


£ £ Go,s\ 1 - rfE(x 2) Mfi„y^-yZ s S Go,s y A~- v 2 E(x 2 ,M,e„,7j P -y 2 )v 


t gl St Q ____ . Q , 

S S G n ^ri-^E(x P ,M,e r ,y,-yp) £ 2H GW 1 - rE{x v ,M, <*„«„- y p ) v ■ ■ . 

11 — 1 v=\ M-I *"1 


S/Sj Q ^ 

_/ 

1 1 


Q 


... ££} G ( ,^ll-v t E(x u M,e r ,y,-yi)v m . . • £ S 

pc=l v = 1 


. . . li S G n ^il- v i E(r 2 ,M,e,,lu-Jh)v m ... £ £ G n ,s y ! l-rE(x 2 ,MA,y*-Th)i 


sjfi r Q 


( 5 ) 


i/s. Q — 

... £ £ c 0 ,svi- 

*=1 »=1 


n 2 E(x , 


,M,dnTu-y P )v m . . . £ 4 G n ^i-v 2 E(x pl MA,y>-yp)v’* 
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The column matrix {j} contains the unknown 
coefficients and the elements are defined as 


floo 

tfoi 


{ 4 = 


n u . 


( 6 ) 


i Q'nm 


From equation (3) the unknown coefficients a nm 
can he obtained by solving the following equation: 

( 4 = [£]-•{«} (7) 

Since the elements of [K}~ 1 can he calculated and 
since the elements of {B} are known from the 
built-in angle-of-attack distribution, the unknown 
coefficients in (.1/ can he evaluated from equation 
(7). Tlie effects of the presence of a central body 
are accounted for by multiplying the slopes in 
{/f} and the downwash due to the hound leg of 
the vortex in the calculation for [/s 1 ]^ 1 in equation 
(7) by the factor [1 + G t2 /y/)]- (See ref. 4.) 

The coefficients determined from equation (7) 
can now he used to calculate the aerodynamic 
parameters c t and c m (appendix C, ref. 4) as follows: 


4 tan atfTT n r, . 9 . 

Cl — - V 1 L-(Ooo + »?‘Uu2+ 


* V Uom) 

+ (« 10 + 1 ? 2 « 12 + • ■ ■ V m dlm)] ( 8 ) 


and 


sir tan a n 9r/ , „ 

0*= -> '\1 — I7‘[(cti 0 +t? «12+ . . . y m ai m ) 

~ (®20 + ?7 2 O22+ • ■ • V m ^2 /«)] (9) 

where the moments are taken about c/4 instead of 
about c/2 as in reference 4. The parameters of 
equat ions (8) and (9) define the rigid-wing aero- 
dynamic loading and will be used to determine a 
corrected angle-of-attack distribution. 


CORRECTED ANGLE-OF-ATTACK DISTRIBUTION 

The change in angle at each control point, due 
to the aerodynamic loading of equations (8) and 


(9), can now be calculated using the relation 


M=[C]{C} +[/>]{ 7’} (10) 


The elements of the column matrix {a<*} are the 
changes in the angle at each control point and 
are defined as 


Ac*! 
Act 2 


{Act} = ' 


(11) 


I At*;,! 


The structural-deflection influence coefficients 
are the elements of the square matrices [C] and 
[D]. The elements C t1 and D i} are the twist at 
the 7th station due to a unit normal load at the 
j th station and the twist at the 7th station due to 
a unit moment at the jth station, respectively. 
The elements of the column matrices {/^} and { T} 
are the integrated loads and moments, respec- 
tively, over the j th spanwise segment; that is, 


i 



A 


and 

L 

T J- V s ' J ] m c 2 dv ' (13) 


where N is the number of segments. The values 
of Ci and c m to be used in equations (12) and 
(13) are obtained from equations (8) and (9), 
respect iv 4v. 

The results of equation (10) (i.e., the changes in 
slopes at the various control points due to the 
aerodynamic loading) arc used to correct the slopes 
of the mean camber line previously used as ele- 
ments of the column matrix {B} in equation (7) 
(before the body factor is applied). These cor- 
rected slopes define the new angle-of-attack distri- 
bution c( erected for aeroelastic distortions. 

ITERATION 

The new angle-of-attack distribution is now 
substituted in equation (7) to determine a new set 
of coefficients a nm . These new coefficients are in 
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turn used in equations (8) and (9) to determine a 
new aerodynamic loading. This loading is then 
used in equation (10) to determine further changes 
in the slopes at each of the control points. The 
process is repeated until there are no further 
changes in the loadings of equations (8) and (9). 
The loading obtained from this iteration process 
is the loading on the deformed wing. Also, the 
angle-of-attack distribution used to obtain this 
loading defines the shape of the deformed surface. 

CALCULATIONS 

Calculations have been made for the wing-body 
combination shown in figure 1. The wing has a 
swcepback of 45° of the 0.25c line, an aspect ratio 
of 4.0, and a taper ratio of 0.15. The wing has an 
NACA 65A206, a = 0 airfoil section at the root and 
varies linearly in thickness to an NACA 65A203, 
a — 0.8 (modified) airfoil section at the 0.50 semi- 
span station; the thickness ratio remains constant 
to the tip. 

This particular wing was chosen because the 
results of references 5 and 6 show that it under- 
went rather large aeroelastic distortions at the 
test conditions. For example, at a Mach number 
of 0.98 and an angle of attack of 4°, the tip twist 



Aspect ratio 4 

Taper ratio 0. 15 

Dihedral, deg 0 

Figure 1 . — Details of wing-body combination. All 
dimensions are in inches unless otherwise noted. 


was —3.54° during tests at 1 atmosphere stagna- 
tion pressure (table VII, ref. 5). In other words, 
the wing tip was operating at only 0.40° when the 
wing-body center line was at 4°. Large aeroelastic 
distortions have a large effect on the spanwise load 
distributions, as noted in reference 6. 

The central body shown in figure 1 was repre- 
sented in the calculations as an infinite cylinder 
with a diameter equal to the maximum diameter 
of the fuselage (3.212 inches). In making the cal- 
culations, the vortex pattern with 12 control 
points and four load lines located as shown in 
figure 2 was used. This vortex pattern divides 
the wing into intervals of 0.1 semispan except in 
the region of the tip where the corrector vortex 
was placed and near the central body where frac- 
tional vortices were employed. In the calculations 
for the downwash factor E, the relations of the 
appendix were used. These relations differ from 
those given in reference 4 in that the bound leg of 
the horseshoe vortices is considered to be swept- 
back to coincide with the swcepback of the load 
lines. As previously mentioned, the slopes used 
in the column matrix {/?} of equation (7) are the 
slopes of the mean camber line at each control 
point. Influence coefficients were obtained for 
this wing from a static calibration and arc pre- 
sented in table I. The influence coefficients were 
obtained along the c/4, c/2, and 3c/4 lines so that 
the changes in slope at all the control points, 
which were also located along these constant- 
chord lines, could he calculated. 

The calculations were programed for the 
IBM 704 electronic data processing machine at the 
Langley Research Center. Generally, five itera- 
tions were necessary at each test condition for the 
aerodynamic and structural forces to be in equi- 
librium. The machine time required was about 
% minute for each test condition. The calculations 
were made at four Mach numbers, A/— 0.80, 0.90, 
0.94, and 0.98, for a stagnation pressure of 1 at- 
mosphere. The results are linear with angle of 
attack so that only enough calculations to define 
the rate of change of the significant parameters 
with angle of attack were necessary. The results 
arc presented at constant angles of attack for con- 
venience in comparing the calculated data with 
the experimental data. 
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COMPARISON OF EXPERIMENTAL AND 
THEORETICAL RESULTS 

Experimental data obtained on the configura- 
tion shown in figure 1 were presented in references 
5 and G and are used herein for comparisons with 
theoretical results. All the experimental data 
used herein were obtained at a stagnation pressure 
of 1 atmosphere. 

SPAN LOAD CHARACTERISTICS 

A comparison of the experimental and theoret- 
ical span load distributions is presented in figure 3. 
At a Mach number of 0.80 (fig. 3(a)) and at angles 
of attack of 0°, 2°, and 4° the theoretical and 
experimental results art* in good agreement. At 
angles of attack of 8°, — 2°, and —4°, the agree- 
ment is poor. Examination of the experimental 
ehordwise pressure distributions of reference 5 
indicates that at these latter angles of attack, a 


leading-edge separation vortex is formed on the 
wing which has significantly altered the character 
of the flow. Viscous effects were not taken into 
account in the theoretical calculations, so predic- 
tion of this type of distribution should not be 
expected. With increases in Mach number to 0.98 
(fig. 3(d)), the agreement at 0° angle of attack 
remains fairly good. However, the agreement 
becomes poor with increases in angle of attack to 
2° and 4°. At the Mach numbers of 0.90, 0.94, 
and 0.98, tin* experimental data of reference 5 
show that the flow over a major portion of the 
upper surface of the wing is supersonic, the extent 
being dependent on both angle of attack and Mach 
number. This flow characteristic is one factor in 
accounting for tin* poor predictions at these Mach 
numbers and angles. If was pointed out in refer- 
ence 7 that the method for performing the 



(a) 47=0.80. (c) 47=0.94. 

(b) 47=0.90. (cl) 47=0.98. 


Fiouke 3. — Comparison of experimental and theoretical spun load distributions for flexible wing. 



<s 


TECHNICAL HE POUT It -58 NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


chordwise integrations by l ho use of replacement for this wing. Figure 5 presents the comparison 

loads (given in ref. 4) breaks down for this wing of experimental and theoretical twist distributions 

plan form as the Mach number approaches 1.0. along the c/4 line at the three higher Mach num- 

This factor may also account for the poor predictions hers. The agreement is still good at angles of 

at Mach numbers close to 1.0. attack from ^2° to 4° at Mach numbers of 0.90 

twist cHAiiACTKitisTics and 0.94 j however, at Al— 0.98, the agreement is 

A comparison of the experimental and theoret- poor at all angles other than 0°. 
ical twist distributions along the c/4, c/2, and 3c/4 < hori wisk centkh-of-prkssukk characterises 

liiu's is presented in ligure 4 for a Mach number of A comparison of the experimental and theoret- 

0.80. In general, the agreement is good at angles ical ehcrdwise center-of -pressure location at an 

of attack from —2° to 4° in each cast*. Figure 4 angle of attack of 4° is shown in ligure 0. In 

also shows that there is little difference in the general, the agreement is good over the inboard 

twist characterist ics at the three ehordwise stations 0.8 of the span up to a Mach number of about 0.90. 



(a) Twist along c/4 line. (b) Twist along c/2 line (c) Twist along 3c/ 4 line. 

Figure 4. — Comparison of experimental and theoretical twist distributions. A/ = 0.80. 
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(a) A/ = 0.90. (b) M =0.94. (c) A/=0.98. 

Figure 5. — ‘Comparison of experimental and theoretical twist distributions along c/4 line, 
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Figure 0. -Comparison of experimental and theoretical 
chordwise center-of-pressure location, a = 4°. 

At a Mach number of 0.80, tlie chordwise pressure 
dist ributions in reference 5 indicate that separation 
was occurring at the 0.95 semispan station at 
a = 4°. This separation accounts for the dis- 
crepancy between the experimental and theoretical 
results at this condition. At a Mach number of 

0.98, the experimental location was rearward of 
the theoretical location over the entire span. 1 he 
distributions of reference 5 show that at this test 
condition, a supersonic distribution developed over 
the trading-edge region and was responsible for the 
rearward center-of-pressure movement. I his rear- 


ward hxation is also responsible for the large dis- 
agreement between the experimental and theoret- 
ical twisr t distributions at ,t/=0.98 in figure 5. 

PRESSURE DISTRIBUTIONS 

Several chordwise pressure distributions have 
been calculated to determine if the good agreement, 
obtained for the other parameters at 4/= 0.80 was 
indicative 1 of good agreement in the pressure dis- 
tributions. Figure 7 presents a comparison ol 
experimental and theoretical chordwise lifting 
pressure distribution at lour spanwise stations and 
at three angles of attack. The* agreement is 
generally good. 

CONCLUSIONS 

An iteration method is presented by which the 
detailed aerodynamic loading and twist charac- 
teristics of a flexible wing with known elastic 
properties may be calculated. The method is 
applicable at Mach numbers approaching 1.0 as 
well as at subsonic Mach numbers. Calculations 
were n ade for a combination of a body and a 
cambered wing with 45° sweepback, an aspect 
ratio of 4, a taper ratio of 0.15, and a thickened 
root section. Comparisons of experimental and 
ealcula ed results at Mach numbers from 0.80 to 
0.98 indicate the following conclusions: 

1 . A a Mach number of 0.80, the calculated and 
experimental loadings and the wing twist distri- 
bution? were in good agreement from 0° to 4° 
angle of attack. At higher positive angles and at 
ncgath e angles, the experimental data indicated 
that se mration was prominent on the wing and the 
agreement between calculated and experimental 
results was poor. 

2. With increases in Mach number, the agree- 
ment between calculated and experimental results 
at 2° and 4° angle of attack became sueceedingly 
poorer. At 0° angle of attack the agreement was 
still good at a Mach number of 0.98. The experi- 
mental data indicated that the poorer agreement 
at angl ‘S of attack of 2° and 4° for the higher Mach 
numbe* was the result of large regions of mixed 
flow which developed on the wing. 

Langley Research Center, 

Nath nai, Aeronautics and Space Administration, 
Langley Field, Va., August t/ f , iftoD. 






APPENDIX 


RELATIONS USED TO CALCULATE DOWNWASH 

In the relations given in appendix A of reference 4 the bound leg of the horseshoe vortices represent- 
ing the wing-body combination is considered to be normal to the five stream. In order to make the 
expression for the dowmvash exact, the sweepback of the wing must be included in deriving the expres- 
sions. Recent unpublished work by John L. Crigler of the Langley Research Center has shown that 
better agreement with experimental results is obtained when the dowmvash is calculated with the bound 
log of the horseshoe vortices swept back to coincide with the load lines. Therefore, these exact down- 
wash expressions were used in the calculations made herein and are presented below. Since a correction 
is applied to the dowmvash associated with the bound leg of the vortex, it is desirable to calculate the 
contribution due to the trailing and bound legs separately. These relations are nondimensional with all 
distances referred to The sum of the contributions of the bound and trailing legs of the vortex 
(width 2.s* P ) at point (/*,//&) to the dowmvash factor at point (s a ^j a ) is given by 

r , 1 1 |~(*« — — ?/«+!)" 

L ’ (x a -7 b ) + t(y b -y a ) 

1 1_ r (x a — ~ ? /a)~hy (x g — y a l) 2 

{ j/b ya I) L (Xq Xft) £ ( //& }fa) 

The contribution of only the trailing legs of the vortex at point (x A ,y & ) and of width 2* c to the dowmvash 
factor at point (x„,?/ a ) is given by 



r t 1_ r (Xg — Xft— 0~h\'(x g — Xft — t) 1 - f fi 2 (yb — Va~\~ i ) ”1 

(Vb Va~\~ 1) L y(Xq — J'b — t) 2j r0 2 {yb — 2/ad“l)" -1 

1 J_ I” (x g ~ y a — l)*n (A 9 ) 

(lU 2 U 1)*.L y(x a — x b -{-t) 2 -\-fi 2 (y b — y a — l) 2 J 

and 


Qba — (fba. k {fba) j b< 


(A3) 


The contribution to the dowmvash factor at the point (x„,?/ a ) ( ue to the vortex at point (x 6) — ?/ & ) for a wing 
with symmetrical loading Al(x b ,y b ) = A£(x&,— y b ) is given by 

r 4- / _ 1 1 F (^~^)+^^~^)+A ; (xg-x & - - 0 2 +d 2 (yft+^ fl +i) 2 n 

Jb. Uba {jy t +y a + J) [_ (?„—**) + < (y*+ j/«) J 

I I” (x g — x ft )-f t (y > — y g )+V(x a — 'X b -\-t) 2 -\-fi 2 (y b -\-y a — 1) 2 H 

(Vb+Va— 1)«.L (Xg —x b ) + t(yt,+ y a ) J 
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and 


f b j _ 1 0 2 H-/? 2 (y ft+y tt ~+ 1 )“~| 

(2/6+2/a 4 -!) XcL \(x a — -Co— O+PWt + Va+iy- 


J 


_ ! 1 Q „— ^+ 0+v(j„--A+0 :t +i9 :g (W i>+.>/ B — l)~ ~| , . r . 

U/>+y.-i)*rL vX'^-^+ 0 2 +^('C+^T)* " J 


Therefore, 


ifba ifba ~T{Jba) f ha (At)) 

In calculating the downwash velocities from the ima^e vortices, it was found that consideration of the 
sweeplmck of the bound log of the vortex had a negligible effect on the results because the width of the 
image vortices is small compared with the width of the wing vortices. Therefore, the slope of the bound 
log was not included in the image calculations and the remainder of the downwash factor was calculated 
by the method of reference 4. 
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TABLE I.— INFLUENCE COEFFICIENTS 


(a) Along c/4 lino 



Wing-twist influence 

coedic'ems due to normal loads at 

Twist 


c/4 

, C,i, dr;, lb, at - — 


measurement 






station, rj 

,' = 0.185 

rj' 0.348 

rj' 0.5)5 

n 0.795 

-0.948 

0. 25 

-0. 0001 

-0. 0002 

0. 0005 

0. 0014 

-0. 0040 

. 40 

0 

. 0002 

- . 001 1 

- . 0032 

0088 

. 00 

1 . 0003 

0002 

0014 

0129 

0215 

. so 

. 0003 

0001 

0015 

- . 0182 

. 0038 

. 05 

. 0005 

0 

— . 0014 

-. 0173 

-. 0950 


Twist 

Wing- twist influence 
r/4. 

eoeOicieuts due to moments 
der/iu-lb, at - 

about 

in visurement 






station, , 

T 0.185 

rj' 0.348 

rj' 0.505 

7)' 0.795 

7)' 0.948 

0. 25 

(). 0001 

0. 0001 

0. 000 l 

0. 0004 

— 0. 0009 

. 40 

. 0001 

. 0004 

. 000 1 

. 0015 

. 0029 

. 00 

. 0002 

. 0000 

. 002 ) 

. 0048 

. 0098 

. 80 

. 0003 

. 0007 

. 003 4 

. 0180 

. 0334 

. 95 

. 0003 

. 0007 

. 004 1 

. 0237 

. 1 130 



(b) Along c/2 line 



Twist 

Wing-twist influence 
c! 4 

coefficien s due to normal loads at, 
, C ih deg lb, at - 

measurement 






station, rj 

V -0.185 

0.348 

rj' 0.505 

V 0.795 

V 0.948 

0. 25 

0 

-•(). 0002 

0. 00 )7 

-0. 0012 

0. 0013 

. 40 

0 

. 0001 

-. 0010 

. 0040 

. 0000 

. 00 

0 

. 0001 

. 0011 

. 0125 

0239 

. SO 

0 

- 0001 

. 00)9 

- . 0107 

. 07 1 0 

. 95 

0 

— . 0001 

00 19 

— . 0150 

. 0835 


Wing-twist influence eoe/liciei ts due 1o moments about 
Twist c/4. />,/, deg/in-lb, at 


• station, , 

V 0.185 

0.348 

7j' -0.505 

,'-■0.795 

V 0.948 

0. 25 

0. 0002 

0. 0002 

0. 0002 

0. 0002 

0 

. 40 

. (1002 

. 0005 

. OOOS 

. 0008 

0 

. 00 

. 0002 

. 0007 

. 0040 

. 0050 

. 0047 

. 80 

. 0002 

. 0007 

. 0042 

. 0233 

. 0383 

. 95 

. 0002 

. 0007 

. 004 1 

. 0258 

. 1453 
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TABLE I. -INFLUENCE COEFFICIENTS —Concluded 
(c) Along 3c/4 lint* 


Twist 

measurement 
station, rj 

Wing-twist influence coefficients due to normal loads at 
c/4, C,„ deg/lb, at — 

7j f — 0.1 85 

t?' -0.348 

V = 0.565 

r,' —0.795 

r,'-- 0.D4H 

0. 25 

0 

-0. 0001 

- 0. 0008 

-0. 0019 

— 0. 0027 

. 40 

0 

0001 

0012 

0048 

0066 

. GO 

0 

— . 0001 

. 0011 

0154 

— . 0305 

. SO 

0 

0001 

— , 0009 

01G5 

0710 

. 95 

0 

0001 

-. 0009 

-. 0153 

— . 0725 


Wing-twist influence coefficients due to moments about 


Twist 

measurement 


c! 4, 

I) a, deg/in -lb, 

at - 


station, rj 

y 0.185 

y = 0.348 

y — 0.5G5 

y 0.795 

V 0.948 

0. 25 

0. 0002 

0. 000 1 

0. 0005 

0. 0004 

0 

. 40 

. 0002 

. 0007 

. 0015 

. 0011) 

. 0023 

. 00 

. 0002 

. 0007 

. 0038 

. 0088 

. 0081 

. 80 

. 0002 

. 0007 

. 0043 

. 0252 

. 0037 

. 95 

. 0001 

. 000G 

. 0050 

. 0252 

. 1434 
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